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Abstract 

Rajagopal and Wilczek have proposed that relativistic nuclear colhsions can gen- 
erate domains in which the chiral condensate is disoriented. If sufficiently large (i.e. 
nucleus sized), such domains can yield measurable fluctuations in the number of neu- 
tral and charged pions. However, by numerical simulation of the zero-temperature 
two-flavor linear sigma model, we find that domains are essentially pion sized. Never- 
theless, we show that large domains can occur if the effective mesons masses are much 
lighter. 



In high multiphcity hadronic colhsions, it is natural to ask if some of the observed pions 
are produced by the coherent decay of a (semi-) classical pion field |l|, 0. Bjorken, Kowalski 
and Taylor have speculated that events can occur in which the pion field is oriented along 
a single direction in isospin space throughout a large fraction of the collision volume 0, 
^. They refer to such domains as "Disoriented Chiral Condensates" (DCC) — disoriented 
because the pion field vanishes in the true vacuum. DCC formation can produce a spectacular 
event structure with some regions of the detector dominated by charged pions and others 
by neutral pions. This behavior may have been observed in Centauro events in cosmic ray 
collisions 0. 

Rajagopal and Wilczek (R&W) have proposed that the nonequilibrium dynamics of 
the chiral transition in relativistic heavy ion collisions can yield large DCC domains. They 
use a linear sigma model to describe the collective chiral behavior, and make the drastic 
assumption that the expansion of the collision fragments "quenches" the field. The initial 
state is taken to be chirally symmetric, as appropriate at high temperature, but its time 
evolution follows the classical equations of motions at zero temperature. R&W argue that 
this quench can amplify the long wavelength pion modes leading to large DCC domains. 

In this Letter we extend the analysis of R&W to estimate the domain size and the 
experimentally relevant fraction of neutral pions TZ^ = n^^o / (n,r+ + ^j^o + n^^- ) . We take all of 
the assumptions of R&W for granted; an analysis incorporating a more realistic description of 
heavy ion dynamics will be presented elsewhere . While our results, like those of R&W, are 
the product of numerical simulations, much of our understanding follows from the analysis 
of Boyanovsky et al. ||^, who study domain growth in the context of inflationary cosmology. 

To identify a nuclear collision in which pion production is dominated by domain formation 
in experiments, one can compare the measured spectra of neutral and charged pions on an 
event-by-event basis ||^ . The probability for a domain to yield a particular fraction of neutral 
pions IZ^ is then P(JZ^) = 1/ (2v/7^3)|p.0||, provided that all isospin orientations are equally 
likely 0, H, |], |Tl|. In contrast, typical hadronic collisions produce a binomial distribution of 
7^3 peaked at the isospin symmetric value of 1/3 (modulo effects related to the isospin of the 
target and projectile). These distributions are very different, particularly at small TZ^. The 
size of domains is of central experimental interest. Nuclear collisions dominated by a single, 
large domain can produce a domain-like PilZ^). On the other hand, events that produce 
many small, randomly-oriented domains result in a binomial distribution. 

To describe the evolution of the classical pion field, we use a linear sigma model that 
characterizes the chiral order parameter of QCD as an 0(4) vector field $ = (cr, tt), where a 
is an isosinglet = 0^ field. The field interacts through a potential A(7f^ + cr^ — f ^)^/4 — iJa; 
at zero temperature the vacuum state is ($) = (/7r,0)- The parameters of the model are 
the coupling A, an external field, H, which breaks the 0(4) symmetry and the vacuum 
expectation value of the sigma field in the absence of of the external field, v. Our primary 
interest is in the physically relevant case of strong coupling, as studied by R&W. One takes 
X = 20, V = 87.4 MeV, and H = (119 MeV)^, so that the pion decay constant is = 
92.5 MeV and the meson masses are 771^^ = 135 MeV and rricr = 600 MeV. To establish the 
role of the ma in domain formation, we compare these results to an extremely weakly coupled 
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case of A = 10 ^, f = 87.4 MeV, and H = (2.03 MeV)^. The pion and sigma fields are very 
liglit in this regime, = 0.3 MeV and m„ = 1.8 MeV, but f^^ retains the measured value. 

Approximating the evolution of the system by a quench, one assumes that the field is 
initially chirally symmetric with ($) = 0, and follows its time evolution using the zero- 
temperature classical equations of motion. To simulate the effect of fluctuations in the 
high temperature initial state, we distribute the fields as gaussian random variables with 
($) = 0, ($2) = v^/A and ($2) = v^/lW following R&W§. Pion domains can form 
in a quench because the chirally symmetric initial state, ($) = (0,0), is unstable against 
small fluctuations in the T = potential @, . In essence, the system 'rolls down' from the 
unstable local maximum of V{^) towards the nearly stable values = /,r^- This process 
is known in condensed matter physics as spinoidal decomposition ^j. Long-lived DCC field 
configurations with jf ^ can develop during the roll-down period. The field will eventually 
settle into stable oscillations about the unique H ^ vacuum, ($) = (/7r,0). Oscillations 
will then continue until interactions eventually damp the motion. In the heavy-ion system 
a domain can radiate pions preferentially according to its isospin content. 

To be more concrete, we write the equations of motion for the Fourier components of the 
pion field in the Hartree approximation |^ : 

= (A(.^-($^))-A:^)vr-i. (1) 

Field configurations with ($) = 0, (<I>^) < v"^ and momentum k < \fXv are unstable in 
that they grow exponentially; modes with higher momentum do not grow. The = mode 
grows the fastest, with a time scale 

Tsp = {X{v^ - ($2))}-^/^ ~ y2/m., (2) 

of order 0.5 fm for A = 20. The exponential growth of the unstable modes continues until 
($^) reaches f ^, when it begins to oscillate about the stable vacuum. Rajagopal and Wilczek 
found that the power oc 7r"^7r"_^ in the low momentum pion modes indeed grows when the 
exact classical equations of motion are integrated for A = 20. What is not clear, however, is 
whether this observed amplification of power is sufficient to create large domains. 

To get a qualitative idea of the domain size, we follow Boyanovsky et al. ^ and estimate 
the contribution of the unstable modes to the pion correlation function. For small times, we 
treat ($2) = v'^/A as a constant, so that 7r^(t) = i'^e*/""'=+%e-*/^% with t^^ = ^J^Xv'^/A + k'^. 
The averages {{AkAk')), {{A^Bk')) over the initial conditions can be computed explicitly 
For t ^ Tsp, we find 

{{TT{x,t)Ti{it))) ^ ((7r(0"',t)'))exp{-a;V8rspt} . (3) 

suggesting that domains grow roughly as S,d ~ y/^Tspt- We therefore expect typical domains 
to reach a size ~ v^Tgp, since they can only grow until the time t ~ Tgp when the field reaches 
the minimum of the potential. If the coupling were weak and the sigma meson light, it would 
take a long time to roll down to the bottom of the potential, so that domains would have 
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plenty of time to grow. However, for the more realistic strongly-coupled case, the roUdown 
is very rapid and the domains are small, perhaps ~ 1.4 fm, which is comparable to the 
Compton wavelength of the pion. 

We confirm this intuition through numerical simulations. Specifically, we simulate the 
model on a three dimensional 10^ x 40 lattice with a spacing a = 1 fm, computing its evolu- 
tion using a simple leap-frog algorithm. To extract domain-size information more easily, we 
employed an asymmetric lattice and studied the pion field averaged over the short dimen- 
sions, TfL{z,t) = J2x,y^{^^yy^yt)/{10)'^. Our simulations qualitatively reproduce the time 
evolution of the momentum space quantity 7i"^7r" ^ reported by R&W for a symmetric lattice 
at A = 20. In configuration space, ttl starts out zero with small gaussian fluctuations. Fig. 
1 shows the profile of vfi at a time t = 100 fm in weak coupling, A = 10~^, and the same 
profile at a time t = 30 fm in strong coupling, A = 20. For weak coupling, regions in which 
the field is slowly varying about some nonzero value — domains — are evident. On the 
other hand, in strong coupling the pion field is oscillating with small amplitude about zero. 
We have also measured the pion correlation function: it is long ranged in weak coupling, but 
drops to zero within a few lattice spacings in strong coupling [0. We also varied the time 
at which correlations were measured, and were never able to find large domains in strong 
coupling. 

We also studied the distribution of the fraction of neutral pions 71^ = ((tt^)^)/ 
Histograms of TZ^ were obtained by evolving 500 independent configurations in time to 
t = 150 and t = 30 fm at weak and strong coupling respectively. These are shown in Figure 
2. In weak coupling the distribution is far from binomial and favors small TZs like the DCC 
distribution Pijls) = 1/(2-^/^). In strong couphng, however, the distribution is clearly 
binomial, peaked about the expected value of 1/3. 

In summary, we find that distributions of the fraction of neutral pions do not reflect 
domain formation for the two-flavor linear sigma model assuming the R&W quench. The 
7^3 distribution can reflect domain formation only when the domain size C,d exceeds the 
system size R. We find that DCC-like behavior dissappers for couplings exceeding A ~ 10~^ 
in our 10^ x 40 lattice, in accord with (3) However, a system the size of a gold nucleus 
can exhibit domain-like behavior for A < 8/(fi?)^ ~ 1. 

A necessary condition for large domains of DCC's to arise is that there is a light particle 
about — the dynamics of a quench does not give us a large distance scale for free. However, 
heavy ion collisions can create a high temperature state that cools over time scales ^ Tgp ~ 
0.5 fm. The quench might therefore not be a a very good approximation, //the equilibrium 
thermodynamics of the QCD phase transition is near a critical point, then near the transition 
temperature there is automatically a light field about, providing a natural mechanism for 
the growth of large domains. In forthcoming papers we will show how such a critical point 
is present in the QCD phase diagram, and analyze the phenomenology of the DCC domains 
produced thereby 

We thank K. Rajagopal for discussions, and for sending us a copy of his thesis. This work 
was supported by a DOE grant at Brookhaven National Laboratory (DE-AC02-76CH00016). 
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Figure Caption 



Figure 1: Snapshot of the pion field, ttl on a 10^ x 40 lattice taken in weak coupling 
(A = 10~^) a,t t — 100 fm and at t = 30 fm at strong (A = 20) coupling. Note how at strong 
coupling the three components of ttl fluctuate around zero. The labeled curves show ttl at 
weak coupling. All three components are non-zero and essentially constant throughout the 
length of the box. 

Figure 2: The distribution of the fraction of neutral pions, TZ^ on a 10^ x 40 lattice in weak 

coupling (A = 10~^) at t = 100 fm and at t = 30 fm at strong (A = 20) coupling. In each case 
500 different initial configurations were used to generate the histograms. At strong coupling 
the distribution is narrowly peaked around 7^3 = 1/3 wheras at weak coupling it is broad 
and clearly giving more weight to small values of TZ^. 
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